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The third and fourth degree collisional moments for d-dimensional inelastic Maxwell models are 
exactly evaluated in terms of the velocity moments, with explicit expressions for the associated 
£Nj ' eigenvalues and cross coefficients as functions of the coefficient of normal restitution. The results are 

applied to the analysis of the time evolution of the moments (scaled with the thermal speed) in the 
free cooling problem. It is observed that the characteristic relaxation time toward the homogeneous 
cooling state decreases as the anisotropy of the corresponding moment increases. In particular, in 
, contrast to what happens in the one-dimensional case, all the anisotropic moments of degree equal 

. to or less than four vanish in the homogeneous cooling state for d > 2. 
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I. INTRODUCTION 



A realistic model capturing the influence of dissipation on the dynamic properties of granular systems consists of 
a gas of inelastic hard spheres (IHS) with a constant coefficient of normal restitution a < 1 For sufficiently low 
densities, the Boltzmann equation for IHS provides the adequate framework to describe the time evolution of the 
one-particle velocity distribution function /(r, v; t) [2J. However, the intricacy of the Boltzmann collision operator for 
IHS makes it difficult to obtain exact results. For instance, the fourth cumulant ai of the velocity distribution in the 
so-called homogeneous cooling state (HCS) is not exactly known, although good estimates of it have been proposed 



[3, 0, For inhomogeneous situations, explicit expressions for the Navier-Stokes (NS) transport coefficients are 



I 

q ■ approximately obtained by considering the leading terms in a Sonine polynomial expansion 

As in the elastic case, part of the above difficulties can be overcome by considering the so-called Maxwell models, 
i.e., models for which the collision rate is independent of the relative velocity of the two colliding particles. Inelastic 
Maxwell models (IMM) have attracted the attention of physicists and mathematicians since the be ginning of the 

century [H, E3, M M M EE El, El El, M M SO, M MMMMM M M M SO, SI M Hj. The 

structure of the Boltzmann collision operator for IMM has the advantage of allowing for the derivation of a number 
of exact properties, such as the high- velocity tails [3, EE, Dj| E3, El, EH and the cumulants [H, E3, H3, 11E IIS Ull in 
r***" ■ homogeneous situations, the NS transport coefficients [27], [3l|, and the rheology under simple shear flow [13|,[2!|- As a 
consequence, it is possible to explore the influence of inelasticity on the dynamic properties in a clean way, without the 
need of introducing additional, and sometimes uncontrolled, approximations. Apart from their academic interest, it 
turns out that the IMM reliably describe the properties of IHS in some situations, as happens in the simple shear flow 
problem and for the NS transport coefficients associated with the mass flux [3l| . Furthermore, it is interesting to 
remark that recent experiments [35[ for magnetic grains with dipolar interactions are well described by IMM. 

The aim of this paper is to contribute to the advancement in the knowledge of exact properties of IMM by evaluating 
1 all the third and fourth degree moments of the Boltzmann collision operator for an arbitrary number of dimensions 
d. The knowledge of those collisional moments, along with that of the second degree collisional moments [27], HI] , 
opens up a number of interesting applications. For instance, one can investigate the temporal relaxation toward the 
HCS, starting from arbitrary initial conditions (not necessarily isotropic), as measured by the lowest degree moments 
(namely, the fourth degree moments) which signal the non-Gaussian character of the asymptotic velocity distribution 
function. This issue will be covered in this paper. 

This paper is organized as follows. In Section |TT] the Boltzmann equation for IMM is presented. Next, the Ikenberry 
polynomials [36| 5^ r u 1 i 2 ...i 1> (V) of degree k = 2r + s are introduced and their associated collisional moments J2r\i 1 i 2 ...i B 
for k = 3 and 4 are evaluated, the technicalities being relegated to an Appendix. The results are applied to the 
relaxation problem of the (scaled) moments toward their asymptotic values in the HCS in Section IIII1 The paper is 
closed in Section ITVl with a brief discussion of the results obtained here. 
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II. COLLISIONAL MOMENTS FOR IMM 

In the absence of external forces, the inelastic Boltzmann equation for a granular gas reads [l[ 

(£> t +vV)/(r,v;i) = J[v|/,/], (2.1) 

where J[v|/, /] is the Boltzmann collision operator. The form of the operator J for IMM can be obtained from the 
form for IHS by replacing the IHS collision rate (which is proportional to the relative velocity of the two colliding 
particles) by an effective velocity-independent collision rate. With this simplification, the form of J becomes [28[ 

J [vi|/, f] = ^-Jdv 2 J da [q-VW/K) - /( Vl )/(v 2 )] . (2.2) 



Here. 



dv/(v) (2.3) 



is the number density, v is an effective collision frequency (to be chosen later), fid — 2n d / 2 /T(d/2) is the total solid 
angle in d dimensions, and a < 1 refers to the constant coefficient of restitution. In addition, the primes on the 
velocities denote the initial values {v'^v^} that lead to {vi, V2} following a binary collision: 

vi = Vl - 1 (1 + a- 1 ) (a ■ E )a, v 2 = v 2 + - (l + a' 1 ) {a • g)5 , (2.4) 

where g = V! — v 2 is the relative velocity of the colliding pair and <x is a unit vector directed along the centers of the 
two colliding spheres. The collision frequency v can be seen as a free parameter in the model. Its dependence on the 
coefficient of restitution a can be chosen to optimize the agreement with the results obtained from the Boltzmann 
equation for IHS. In particular, to get the same expression for the cooling rate as the one found for IHS (evaluated in 
the local equilibrium approximation) one takes the choice [271 ] 

d + 2 d-i [T , oc s 

„ = — »o, v«= ^ {d + 2) ™ Vm' (2 - 5) 

where a is the diameter of the spheres. Note that, in any case, the results derived in this paper will be independent 
of the specific choice of vq . 

A useful identity for an arbitrary function h(v) is given by 

J[h] = J dv 1 / l (v 1 )J[v 1 |/, f] = -^-fdv 1 J dv 2 /(vi)/(v 2 ) J da [h(y'{) - ft( Vl )] , (2.6) 

where 

v'/ = v 1 -i(I + a)(^-g)^ (2.7) 

denotes the post-collisional velocity. If h(v) is a polynomial, then 

M[h] = J dv/i(v)/(v) (2.8) 

is its associated velocity moment and J[h] is the corresponding collisional moment. 

In the case of Maxwell models (both elastic and inelastic), it is convenient to introduce the Ikenberry polynomials 
[H[ Y2 r \i 1 i 2 ...i g (y) of degree k = 2r + s, where V = v — u(r) is the peculiar velocity, u(r) being the mean flow velocity 
defined as 

u=~ [ dw/(v). (2.9) 



The Ikenberry polynomials are defined as Y 2r i ili2 _ ia (V) = V 2r Yi li2 ,,,i s (V), where iiii 2 ...i s (V) is obtained by sub- 
tracting from Vi x Vi 2 . . . Vi s that homogeneous symmetric polynomial of degree s in the components of V such as to 
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annul the result of contracting the components of 5^ 1 j 2 ...i s (V) on any pair of indices. The polynomials functions 
^2r\iii 2 — ieflf) of degree smaller than or equal to four are 

y |o(v) = i, %(v) = y i , (2.10) 

^2|o(V) = ^ 2 , y |y(V) = W,-^ 2 % (2.11) 

Y 2 \i{V) = V 2 Vi, Y 0]ijk {V) = ViVjVk ~ j^V 2 {ViSjk + VjS ik + V k 5 tJ ) , (2.12) 

1 



Y 4l0 (V) = V\ Y m (V)=V 2 [V i V j --V 2 S ij ), (2.13) 



Yo\wa<y) = ViVjVkVt - j^V 2 QWjSkt + ViV k 6 je + V t V t 5 ]k + VjV^u + VjV e S ik + V k VA,) 

V A {8ij5 k e + 6 ik Sji + SuSjk) 



(d + 2)(d + A) 

= ViViVM - [Y m (V)5 ki + Y 2llk (V)6 j£ + Y 2lu (V)S jk + Y m (V)S u 

+Y m (V)S ik + F 2 | M (V)%] - — — — V 4 (S i:j 6 k e + S ik S je + S u S jk ) . (2.14) 

Here we will use the notation M 2r \i 1 i 2 ...i s = -M [Y 2r \ ili2 ^ Aa ] an< ^ Jir\i^...i s = J\Y2 r \i 1 i 2 ...i s ] f° r the associated 
moments and collisional moments, respectively. Note that Af | = n, Jq|o — (conservation of mass), M u = (by 
definition of the peculiar velocity), Joij = (conservation of momentum), and Af 2 |o — pd/m, where p — nT is the 
hydrostatic pressure, T being the granular temperature. Moreover, Mquj = (Pij — pSij)/m, where Pij is the pressure 
tensor and M 2 u = 2qi/m, where q is the heat flux vector. The moment M 2 \q, the number density n, and the flow 
velocity u are the hydrodynamic fields, while the moments M$uj and M 2 \i constitute the momentum and energy 
fluxes, respectively. The remaining third degree moments M uj k and the moments of a degree k > 4 are not directly 
related to the hydrodynamic description, but they are useful to provide information about the velocity distribution 
function. In particular, the moment M 4 | is related to the fourth cumulant a 2 as 

o 

171 

a2 = d(d + 2)nT^ - 1 > (2 - 15) 

while the moments Moiyfe, Mqu^, and M 2 uj measure the degree of anisotropy of the velocity distribution. 

As in the elastic case, the mathematical structure of the collision operator (|2.2p implies that a collisional moment 
of degree k can be expressed in terms of velocity moments of a degree less than or equal to k. More specifically, 
the choice of the polynomials Y 2r \ s (V), where we have introduced the short-hand notation s = i\i 2 . . .i s , yields the 
following structural form for the collisional moments J 2r \ s - 

J 2 r\s — ~ v 2r\sM 2 r\s + X r i r n | g'g"gM 2r i | gi M 2r ii | g « , (2-16) 

r' ,r" ,g' ,g" 

where the dagger in the summation denotes the constraints 2(r' + r") + s' + s" = 2r + s, 2r' + s' > 2, and 2r" + s" > 2. 
Since the first term on the right-hand side of Eq. (|2. 16[) is linear, then v 2r \ s represents the eigenvalue of the linearized 
collision operator corresponding to the eigenfunction Y 2r \ s (V). 

Let us now display the explicit expressions for the collisional moments J 2r \i 1 i 2 ...i a for k = 2r + s < 4. We start with 
the second degree moments. 

A. Second degree collisional moments 

The second degree collisional moments were already evaluated in Ref. 1271 . They are given by 

J 2 \o = — v 2 \oM 2 \ , J \ij = — vo\ 2 M \ij, (2-17) 
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where the expressions for the eigenvalues v 2 \a and Vq\ 2 are 

^ 2 |o = ^ (1 - « 2 ) Ho, (2.18) 

(l + a)(rf+l-a) (1 + a) 2 
^o|2 = = 2 I° H 4 °' ^ ' 

The quantity i/ 2 |o is not but the cooling rate, i.e., the rate of change of the granular temperature due to the inelasticity 
of the collisions. The eigenvalue vo\2 is the collision frequency associated with the NS shear viscosity and reduces 
to vq in the elastic limit. The second equality in Eq. (|2.19[) decomposes f |2 into the part inherent to the collisional 
cooling plus the genuine part of the momentum collisional transfer. As shown below, a similar decomposition can be 
carried out for the eigenvalues v 2r \ s . 

B. Third degree collisional moments 

The evaluation of the third degree collisional moments J 2 \i and Jo\ijk is performed in the Appendix. The results 
are 

J%\i = -v 2 \iM 2 \i, Jo\ijk = _ fo|3-Wo|tjk) (2.20) 

where 

(1 + a) [5d + 4-a(d + 8)} 3 (1 + ^(d- 1) 

^ ~ M V ° 2^° + 4d U0 > (2 - 21) 

3 

^o|3 = 2^o|2- (2.22) 

Equation (|2.21j) was first obtained in Ref. [23- The eigenvalue v 2 \i is the collision frequency associated with the NS 
thermal conductivity. It is interesting to note that (1*211 — § ^2|o)/( I/ o|2 — ^2|o) — (A ~ 1)/^ which generalizes the 
simple relationship, holding for elastic Maxwell models, between the collision frequencies associated with the thermal 
conductivity and the shear viscosity. An even simpler extension is provided by Eq. (|2.22[) . 

C. Fourth degree collisional moments 

The fourth degree collisional moments are also worked out in the Appendix. They can be written as 

Ja\o = -^4|oA^4|o + Aira -1 Mf| - X 2 n^ 1 M Q \ lj M \ ji , (2.23) 



(2.24) 



Jolijkl — — ^014-^01 ijM + X5TI 
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M \ijM \ki + M \ ik M y e + M Q \ ie M \ jk - — — {M \ ip Mo\ P j5ke 
p Mo\pe^jk + M \j p Mo\pf.6u + i 
Mo\ P qM \ qp (Sijdke + Sik6j£ + SuSjk) 



+M \ ip M \ pk Sj£ + M \ ip M \ pe Sjk + M Q \j p M \ p )-6u + M \jpM \ pi 5ik + M \k p M \ p z6ij 
2 

+ (d+2)(d + 4) 



(2.25) 



In Eqs. (|2.23I) - (I2.25|) . the usual summation convention over repeated indices is assumed. The collision frequencies (or 
eigenvalues) v 2r \ s and the cross coefficients Xi are given by 

(1 + a) [I2d + 9 - a(4d + 17) + 3a 2 - 3a 3 ] 

16d V ° 

(1 + a) 2 Ud- 7 + 6a -3a 2 ) 
= 2^ 2 | + ^ V 1( , d } -vq, (2.26) 
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FIG. 1: (Color online) Plot of the (reduced) fourth degree eigenvalues i/Jj, and i/^i (top panels) and of the shifted 

eigenvalues W014, ^2|2j an d w 4|o (bottom panels) as functions of the coefficient of restitution. The left and right panels correspond 
to d — 2 and d = 3, respectively. 



(1 + a) [id 2 + 31d + 18 - a(d 2 + Ud + 34) + 3a 2 (d + 2) - 6a 3 ] 
1/212 = 8d(d + 4) ^ 

(1 + a) 2 [3d 2 + 7d - 14+ 3a(d + 4) - 6a 2 ] 

= 2i/2 ' + 8d(dTI) 



(2.27) 



^0|4 



(1 + a) [2d 3 + 21d 2 + 61d + 39 - 3a(d + 3)(d + 5) + 3a 2 (d + 3) - 3a 3 ] 



2^2|C 



2d(d + 4)(d + 6) 
(1 + a) 2 [d 3 + 9d 2 + 17d - 9 + 3a(d + 4) - 3a 2 ] 
2d(d + 4)(d + 6) 



(2.28) 



(1 + a) 2 (d + 2) (4d - 1 - 6a + 3a 2 ) 
Al = 16d 2 V(h 



(2.29) 
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(l + g) 2 (l + 6a-3a 2 ) 

A 2 = -vq, 2.30) 

8a 

(1 + a) 2 [d 2 + bd - 2 - 3a(d + 4) + 6a 2 ] 

A 3 ^-p vo, (2.31) 

(l + a ) 2 [2-rf + 3 a (rf + 4)-6« 2 ] 

A4 4dJd + 4) **' (2 ' 32) 

_ (1 + a) 2 [rf 2 + 7rf+9-3 a (rf + 4) + 3 a 2 ] 
As " 2d(d + 4)(d + 6) "°- (2 - 33) 

Equations 12.26(1 and (12.29|) coincide with the results of Ref. H3. 

We have checked that, in the elastic case (a = 1) and for three-dimensional systems (d = 3), all the expressions 
reported in this Section reduce to known results [3g, In the one-dimensional elastic case {d = 1, a = 1), the 
gas behaves as an ideal gas because a collision is equivalent to exchanging the labels of both colliding particles. As 
a consequence, J2 r \s = 0. It is easy to check that Eqs. (|2.17p - (|2.33p are consistent with this property since the 
coefficients affecting the non- vanishing moments are zero, i.e., v 2 \i = v i\a = Ai = 0. Moreover, for one-dimensional 
inelastic gases (d = 1, a < 1), our expressions for i/ 2 \o, t>2|i, ^410) an( I Ai agree with the results derived by Ben-Nairn 
and Krapivsky fl2l |, who obtained the exact expressions for all the collisional moments, namely J 2r \o an( i Jir\x- 

While the a dependence of the second and third degree eigenvalues, Eqs. (|2.18[) . (|2 . 19[) . (|2.21[) . and (|2 . 22[) . is 
relatively simple, that of the fourth degree eigenvalues (|2.26[) - (|2.28p and the cross coefficients (|2.29[) - (12.33[) is more 
involved. Figure [T] shows the a dependence of the (reduced) eigenvalues Vm , ^212' anc ^ ^014' wnere v 2r \ s = v 2r\s/v{), 
and the shifted eigenvalues w 4 | , oj 2 \ 2 , an( i w o|4i where we have called uj 2r \ s = ^ 2r | s ~~ ( r + s /2)^2|o» d = 2 and d = 3. 
While ^q| 4 decays monotonically as the inelasticity increases, the other two eigenvalues v%\ 2 an d v*l\o start growing, 
reach a maximum, and then decay. The maximum value of u 2 , 2 occurs at a ~ 0.40 for d = 2 and at a ~ 0.67 for e? = 3. 
In the case of z^| , the maximum occurs at a ~ 0.18 and a ~ 0.30 for <i = 2 and o? = 3, respectively. However, when 
the part associated with the cooling rate is subtracted from the bare eigenvalues, the resulting shifted quantities u>mq, 
uj 2 \ 2 , an d w o|4 exhibit a monotonic behavior. As shown in Section [Ml these shifted eigenvalues are the relevant ones 
in the time relaxation of the scaled moments in the free cooling problem. Therefore, the decrease of uj^q, u> 2 i 2 , an( I 
w |4 implies that the characteristic relaxation times of the (scaled) fourth degree moments toward their asymptotic 
values increase with dissipation. 

It is instructive to compare the fourth degree eigenvalues with the second and third degree eigenvalues. In the 
elastic case, one has lo^q = lu 2 \i < ^o|2 — w 2|2 < w o|3 < w o|4 f° r d — 2 and U4\ = u 2 \i < u>o\2 < w 2| 2 < w o|3 < w o|4 
for d = 3. We have observed that inelasticity breaks the degeneracy lj 4 \q = ui 2 \i for both dimensionalities (yielding 
W4|o < u>2\i) and the degeneracy u>q\ 2 = u 2 \ 2 for d = 2 (yielding io 2 \ 2 < ojq\ 2 ). The inelasticity also affects the ordering 
of the eigenvalues: for d = 2 one has cj 4 | < u 2 \i < ^ 2 \ 2 < ^o|2 < w o|3 < ^014 if 0-17 < a < 1 and w 4 | < £^211 < 
(jj 2 | 2 < ^o|2 < w o|4 < ^013 if < a < 0.17; for d = 3 the ordering is lu^ < ct> 2 |i < ^o|2 < ^212 < ^013 < w o|4 if 
0.43 < a < 1 and W4|o < ^2|i < ^2| 2 < ^o|2 < w o|3 < ^014 if < a < 0.43. Since, except w 4 | , these quantities are 
related to moments which vanish in isotropic states, the fact that w 4 |o is the smallest one implies that (as expected 
on physical grounds) the characteristic time needed to achieve an isotropic state is shorter than the one needed to 
reach the asymptotic state. 

Let us consider now the (reduced) cross coefficients A* = Xi/vo (i = 1, . . . , 5), which measure the coupling of the 
second degree moments to the evolution of the fourth degree moments. Their dependence on dissipation is shown in 
Fig. O It is apparent that the effect of inelasticity on A* is more pronounced than on the fourth degree eigenvalues. 
For elastic collisions, A3 = < A5 < A| = A 2 = X\ for d = 2 and A5 < A3 < A4 = A 2 < A* for d = 3. This 
ordering changes with inelasticity. Moreover, A*, A 2 , and A4 significantly decrease with increasing dissipation, A3 has 
a non-monotonic behavior, and A5 is nearly constant. Note that the coefficient A4 does not actually play any role 
in d = 2 since the combination Mqu^Mq^j ~ \^o\kl^oUk^ij appearing in the collisional moment J 2 u, cf. Eq. (12.241) . 
vanishes in the two-dimensional case. 



III. RELAXATION TO THE HOMOGENEOUS COOLING STATE 



The results derived in the preceding Section can be applied to several interesting situations. Here we will consider 
the most basic problem, namely the time evolution of the moments of degree less than or equal to four (both isotropic 
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and anisotropic) in the homogeneous free cooling state In that case, the Boltzmann equation (|2.ip becomes 

dtf(v,t) = J[v\fJ], (3.1) 

which must be complemented with a given initial condition /(v,0). Since the collisions are inelastic, the granular 
temperature T(t) monotonically decays in time and so a steady state does not exist. In the context of IMM, it has been 
proven [HLHl! that, provided that /(v, 0) has a finite moment of some degree higher than two, /(v, t) asymptotically 
tends toward a self-similar solution of the form 

/(v,i) -> n[v (t)]-U(V/v (t)), (3.2) 



where v (t) = sj2T(t)/ m is the thermal speed and 0(c) is an isotropic distribution that is only known in the one- 
dimensional case [13] • According to Eq. (| 3 - 2 [) , the scaled moments 

M* r|s _(t) = n-> (t)r (2r+s) M 2r | s -(i) (3.3) 

must tend asymptotically to 



M 2 * r|s -(i) -> ii2r\s = J dcY 2r | s -(c)0(c). (3.4) 

Due to the isotropy of 0(c), then /i2r\s = unless s = 0. Moreover, it is known that the scaled distribution 0(c) 
exhibits an algebraic high velocity tail [3, EH EH [3 of the form 0(c) ~ c~ d ~' y< - a \ so that the moments /x 2r |o diverge 
if 2r > 7(a). The quantity 7(a) obeys a transcendental equation whose solution is always 7(a) > 4, except for 
d = 1. Consequently, for any value of a and <i > 2, the scaled moment M^ (t) goes to a well defined value /i4|o, 
while the remaining scaled moments of degree equal to or less than four are anisotropic (except of course M^ Q = 1 
and -Mjhq = d/2) and so they tend to zero. The main goal of this Section is to analyze in detail the relaxation of the 
second, third, and fourth degree moments (both isotropic and anisotropic) toward their asymptotic values. 
Taking velocity moments in both sides of Eq. (|3.1|) one has 

dtM 2r \s = J2r\s- (3.5) 

In particular, 

d t M 2 \o = -^ 2 |o^2|o- (3-6) 
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Since M 2 \o — dnT/m, Eq. (]3.6|) is the equation for the time evolution of the granular temperature and V2\o is the 
cooling rate. The solution of Eq. (|3.6[) is 

T(t) = 3 ff , (3.7) 

[1 + i*,|o(0)t/2] 

where T(0) is the initial temperature and z^loCO) oc T 1 / 2 (0) is the initial cooling rate. Equation (|3.7[) is not but Haff's 
law Q. 

Let us consider now the scaled moments Q3.3p . In that case, from Eqs. (13. 5|) and (|3.6[) one simply gets 



where J£, 5 = J 2 , |s/^o«Wo r+S and 



a T M* r|s - = J* r|5 + ^±l^ |0 M* r|g , (3.8) 



r 



= / &t'v {t') (3.9) 



i) 



measures time as the number of (effective) collisions per particle. The effect of the second term on the right-hand 
side of Eq. (|3.8[) is to shift the eigenvalues i/£ . to UJ 2r\s — v 2r\s ~ ( r + s /2) v 2\o- ^or instance, 



d T M* MQ = -u, 4|0 M 4 * |0 + K^- KM*^M* Ui . (3.10) 



A. The one-dimensional case 



In the one-dimensional case, the only scaled moments of degree equal to or less than four are (apart from M*^ = 1 
and M*| = |) MZ and MJ . Their evolution equations are 

d T M* ]x = 0, (3.11) 



d T Mt ]Q = -{l-a 2 Y 



The solution of Eq. (gUg) is 



^4*|0 



(3.12) 



M 4 * |0 (r)=M 4 * |0 (0)eAd- 2 ) 2 - + £ 



3 ^(l- Q 2 ) 2 r _ j 



(3.13) 



This solution shows that the (scaled) fourth degree moment monotonically increases with time, i.e., iamq = oo. This 
is consistent with the exact HCS solution found by Baldassarri et aZ. [T^|. namely 



(c) 



2 3/2 



1 



7T (1 + 2C 2 ) 



(3.14) 



On the other hand, Eq. p. lip shows that M^ x (t) = M^ x (0), i.e., if the initial state is anisotropic with M2| x (0) 
then one has M 2 \ x (t) = M 2 | 2; (0)[T(i)/T(0)] 3 / 2 . The constancy of M^, x implies that any initial anisotropy does not 
vanish in the scaled velocity distribution function for long times. As a consequence, while the distribution (|3. 14|) 
represents the asymptotic form <f>(c) for a wide class of isotropic initial conditions, it cannot be reached, strictly 
speaking, for any anisotropic initial state. Whether or not there exists a generalization of (|3.14p for anisotropic 
states is, to the best of our knowledge, an open problem. Since the symmetry of the distribution (|3.14[) implies that 
(c x ) — but the average (|c| 3 ) diverges, a small correction to the form (|3.14p could accommodate a finite 



M, 



2\x 

value of (c 3 ) 
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B. The two-dimensional case 



As is known, one-dimensional systems are generally not very realistic and so they can exhibit peculiar properties. 
However, two-dimensional systems are usually considered as representative of the features present in real systems. 

For the sake of simplicity, in the remainder of this Section we will consider the two-dimensional case. The set of 
independent moments of second, third, and fourth degree will be taken as 

{M^ ]xx ,M* lxy }, (3.15) 

{M* w M* {y , M* |xxx , M* |xx J, (3.16) 

{M 4 j , M 2 \ xx , M 2 \ xy , M \ xxxx , M Q \ xxxy }, (3-17) 
respectively. The remaining moments are simply related to the above ones as M^ yy = —M^,, MZ yy — — MZ xxx , 

M 0\yyy = "^oW M 2\yy = M 0\yyyy = = M oW M 0\ Xy yy = 'M^y From Eq. ftS), it is 

easy to obtain the time dependence of the (scaled) second and third degree moments: 

M^ ij (r)=M* lij (Q)e-^, (3.18) 

Ml {i (r) = M^Oje-*^, M* m (r) = M* m (0)e'^ . (3.19) 

In the case of the fourth degree moments, one has to deal with inhomogcncous linear differential equations involving 
the second degree moments. The solutions are 



M 4 * |0 (t) = M| |0 (0)e-^+-^(l-e 



w 4|0 

2X% 



v 2 



2Wq|2 — w 4|0 



^oL(O) + M t JO) (e-»^ - e ~ 2 ^) , (3.20) 



M* m {r) = M 2 * to .(0)e-^ + % M * |y (0) ( e "^i^ - e ~^) , (3.21) 

' ' W |2 - W 2 | 2 ' 

MS\ xxxx (t) = M* ]xxxx (0)e-^r + 3 X* 5 [m *! 2 xx (0) - M*? xy (0)} ( e -*l' T - e ~^) , (3.22) 

1 2 2W | 2 - W |4 L 1 1 J 

3A* 

M^ xxxy (r) = M* lxxxy (0)e-^ + 5 M * |xx (0)M o * |xj/ (0) (e"^^ - e ~^ r ) . (3.23) 

^0|2 — W |4 

Equations (|3.18p - (|3.23|l show that all the moments, except MJq, tend to zero for sufficiently long times. The asymp- 
totic expression of M^ Q is 

X{ _ r) 7-6a + 3a 2 

w 4|0 



M 4 * |0 -> A*4|o = — = 2-— r-j, (3.24) 



which agrees with previous results [27] . 

As an illustration, let us analyze the time evolution of the scaled fourth degree moments (|3.17|) for the following 
initial anisotropic distribution: 

77 

/(v, 0) = - [6(Y - Vi) + S(V - V 2 ) + 5(V - V 3 )] , (3.25) 



where Vi = u (0)x, V 2 = (u (0)/i/2)y, and V 3 = -Vi - V 2 . Here, v (0) = y/2T{&)/m is th e initi al thermal speed, 
where the initial temperature T(0) is arbitrary. Figure [3] shows the evolution of the moments (j3.17[) for two values of 
the coefficient of restitution: a = 1 (elastic system) and a = 0.5 (strongly inelastic system). It is quite apparent that 
the number of collisions needed to reach the HCS values increases with the inelasticity, as expected from Fig. [T] In 
the particular case of a = 0.5, the relaxation times are about twice the ones corresponding to a = 1. Moreover, since 
w | 4 > w 2 |2 > w 4|0i we observe that the moments M*^ ke tend to zero more rapidly than the moments M*^, and that 
the isotropic moment M£, Q reaches its asymptotic value more slowly than the anisotropic moments. 
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FIG. 3: (Color online) Time evolution of the scaled fourth degree moments in the free cooling of a two-dimensional IMM with 
a — 0.5 (solid lines). The time evolution to equilibrium of an elastic system (a = 1) is also shown (dashed lines). In both 
cases the initial state is described by Eq. (|3.25|l . The horizontal dotted lines in the bottom panel indicate the corresponding 
asymptotic HCS values flM. 



IV. DISCUSSION 

As Maxwell already realized [38| , scattering models where the collision rate of two particles approaching each other 
with a relative velocity g is independent of the magnitude of g allows one to evaluate exactly the collisional moments 
without the explicit knowledge of the velocity distribution function. In the conventional case of ordinary gases of 
particles colliding elastically, Maxwell models are useful to find non-trivial exact solutions to the Boltzmann equation 
in far from equilibrium situations [37| • Needless to say, the introduction of inelasticity through a constant coefficient 
of normal restitution a < 1 opens up new perspectives for exact results, including the elastic case as a special limit 
(a = 1). This justifies the growing interest in IMM by physicists and mathematicians alike in the past few years. 

The choice of the Ikenberry polynomials [36| Y 2r \ s of degree 2r + s allows one to express the corresponding collisional 
moments J2 r \s m the form (|2. T6[) : an eigenvalue —V2r\s times the velocity moment M 2r \s plus a bilinear combination of 
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moments of degree less than 2r + s. In particular, u 2 \o is the cooling rate of the gas. In this paper we have evaluated all 
the third and fourth degree collisional moments of the IMM defined by Eq. (|2.2[) . In that context, the results are exact 
for arbitrary values of a and apply to any dimensionality d. Known results are recovered for three-dimensional elastic 
systems [36|, H3] and for one-dimensional inelastic systems [l^]. We have observed that some of the eigenvalues v 2r \ s 
do not have a monotonic dependence on a, while the shifted eigenvalues v 2r \ s ~ ( r + s / 2 ) 1 o monotonically decrease 
with increasing inelasticity. Moreover, given a value of a and a degree 2r + s, the eigenvalues ^2r|s increase with s. 
This means that the larger the anisotropy of a moment M 2r \s the higher its collisional rate of change. Although the 
above observations are based on the moments of degree 2r + s < 4, we expect that they extend to moments of higher 
degree. 

As a simple application of the results derived in Section [TlJ we have studied the time evolution of the moments of 
degree equal to or less than four in the free cooling state, in which case the evolution of the moments scaled with the 
thermal speed is governed by the shifted eigenvalues. An interesting feature of the one-dimensional case is that the 
heat flux q x = (m/2)M 2 \ x , when scaled with the thermal speed, does not change in time, so that an initial anisotropic 
distribution cannot evolve toward an asymptotic isotropic distribution. Thus, the exact solution found by Baldassarri 
et al.[l7\ does not play a universal role, at least in a strict sense, unless the initial distribution is isotropic. On the 
other hand, we have found that all the anisotropic moments of degree equal to or less than four vanish in the long time 
limit for d > 2. However, this does not preclude the possibility that anisotropic moments of higher degree diverge for 
a sufficiently small. We plan to explore this possibility in the near future. 

The explicit results provided in this paper can be useful for studying different problems. An important application 
is the exact derivation of the Burnett order constitutive equations for IMM, with explicit expressions of the associated 
transport coefficients as functions of d and a. This is possible because the determination of the Burnett order pressure 
tensor and heat flux requires the previous knowledge of the third and fourth degree collisional moments to Navier- 
Stokes order. Another interesting problem is the so-called simple or uniform shear flow, which is an intrinsically 
non-Newtonian state [39| . Apart from the rheological quantities, the results derived here allows one to analyze the 
time evolution of the fourth degree velocity moments toward their steady state values and investigate their 
possible divergence, in a similar way to the analysis carried out in the elastic case [4l|. Finally, the generalized 
transport coefficients characterizing small perturbations around the simple shear flow have been determined [42| and 
compared with those previously obtained for IHS [43] from a model kinetic equation. 
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APPENDIX: EXPLICIT EVALUATION OF THE COLLISIONAL MOMENTS 



In this Appendix we give the details of the derivation, by using the property (|2.6p . of the collisional moments 
J2r\s = c/[^2r|s] associated with the Ikenberry polynomials of third and fourth degree. To carry out the calculations 
we will need the angular integrals 



J da {a ■ gf r +% = B r+1 g 2r 9l 



(A.l) 



da (a ■ g) 2 ^ = ^.9 2(r - 1} (2rg i9j + , 



(A.2) 



/D 
da {a ■ g) 2r+1 did,jd k = — + [ 2r 9i9j9k + g 2 {kj9k + kk9j + ^jk9i)] , (A.3) 



/ 



4r(r - l)g 2{r 2) gtgjgkge 



{2r + d)[d + 2(r + 1)} 
+g 2{r ~ 1 ">2r (gigjSke + gig^jt + 9i9^jk + 9k9j^ii + 9e9j&ik + 9k9lhj) 
+g 2r (SijSki + S ik S je + SjkSu)] . 



(A.4) 



Here, the coefficients B r are Q 
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13, = I da (a ■ = ^-^ r ^| r +^ , (A.5) 



1. Third degree moments 



VlWjVH = VxiVtjVtk - ^^(a- ■ g) (VuV^dk + V u V lk aj + V^V^at) 

\ 2 / 



We start by noting that the collision rule (|2.7[) implies that 

1 + a 

T 

+ ( "~7~~ ) ' g) 2 (ViiCTjfffc + Vijdidu + Vikdidj) - ( 1 + a ) (ct ■ g) 3 ^^^. 



Next, making use of Eqs. (|A. (]A.3|) . one gets 

J^uVyVifc] - -^^^ i( d + 2){9iVijVu + 9jV u V lk + gkVuVu) 



1 ~t a [2(gigjVik + gigkVij + gjgkVu) + (g 2 (VuSjk + VijS ik + VikSij))] 



where the brackets are defined as 



JWuVtjVik] 



2d 2 



3(d + 1 - a)M \ ijk + 5d+ ^ d "^ + 8 ^ (M 2 \iSjk + M 2Yj 8 lk + M 2{k S K 



(A.6) 



(A.7) 



<MVi,V 2 )) = ^ / dVi J dV 2 / l (V 1 ,V 2 )/(V 1 )/(V 2 ). (A.8) 
and we have taken into account that (gigjg k ) = (g 2 gi) = 0. It is easy to get 

n{giVijVi k ) = n{gigjVik) = M Q \ ijk + (M 2 \iS jk + M 2 \ 3 5 lk + M 2 \ k 5 l3 ) , (A. 9) 

n(g 2 V u } = M 2 \i- (A.10) 

Therefore, 

vq l + a 



(A.11) 



If one makes j = k and sum over j one gets the first equality of Eq. (|2.20p with !/ 2 |i given by Eq. (|2.21[) . Also, by 
subtracting (J 2 \i5j k + J 2 \j5i k + J 2 \k^ij) /(d + 2) from both sides of Eq. (IA.11I) one gets the second equality of Eq. 
(l2~2"0)) with f | 3 given by Eq. (l2~2"2"|) . 

2. Fourth degree moments 

Now the starting point is the collision rule 

VuV{jV"kVu = VuVijVuVu - ^-^(ff • g) (ViiVyVik?/ + VuVijVuGk + V u V lk V ie S 3 + V^V^Vim) 
l + a\ 2 ,^ , 2 



+ ( - ) (S- ■ g) (VltVij-CTfeS 1 / + VuVikdjcrt + VuVudjdk + VijVikVidi + VijVu^k 



11 N 3 

1 + a 



+VifcVl«<7 i (T 7 -) - I — - — ) (<r ■ g) J {Viidjd k di + Vi^dkde + Vxk^iSj^i + Vu^j^k) 
1 + 



+ — r— (cr ■ gydidjd k dz. (A.12) 
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After integrating over 

J[ViiVijV lk Vu] = - 



l+ " 1 (d+2)(d + 4)(d + 6){g i V lj V lk Vir {: '' K 



2d(d + 4)(d + 6) 2 
1 + a 



-(d + 4)(d + 6) 



2(g igj V lk V ie + ■■■) + {g 2 {V li V ij 6 u + ■ • •)> 



+3(d + 6) 



1 + a 



-3 



1 + a 



%{gi9]9kgt) + ^(g 2 (9t9jSkt.+ • (5) -)) + (fl ,4 )(<%fe+ ■■■) 



(A.13) 



W 



where • • ■ denotes the s terms obtained from the canonical one by permutation of indices. Making k — I and summing 
over k we obtain 

J^VuVu) = ~ 2d( 7° \ U {(d + 2)(d + 4)(2(g • Wx)V u V Xj + V 2 (g i V lj + 9j V u )) 

-i±^(d + 4) [2(^.9,9, + 2(g • VxJOhVy + Si V M )) + ((d + 6) fl V M V y + <? 2 V?%>] 



1 + a 



[(d + 6)<<?%14, + ffi^ii)) + 2((g • Vx){2g i9j + g 2 8 t] ))} 



-3[^) (4.g 2 ff ,9 J+ .9%> 

Next, by taking i = j and summing over i, Eq. (|A.14|I yields 

J[vt] = _^l^| 4(d+2)(fe . Vi) ^ > _ 2 l^ [(d + 4) ^ >+4{(g . Vi)a>] 



(A.14) 



-12 



1 + a 



((g-V!) 5 2 )-3 



1 + Q' 



(3 4 ) 



(A.15) 



2 y u ° ±/a ' v 2 

Now we express the averages in terms of the Ikenberry moments. Let us consider first the four-index averages: 

(A.16) 



nigiVxjVikVu) = M \ ljkf + [M 2 \i 3 5 ke + - (5) - J + ^ M 4 | — ) ■ 



d(d + 2)" 



^(ffiffj^ifeVi^) = M Q \ ljM + ^ 4 [M 2 \ij5u+ ) + ^ ; o ^ M 4 | ( 8ijS ke + ( 



d(d + 2)' 



+n l M \ ij M (s \ k t + ^-M 2 | (M |ij^ + M |fc«<%) + -^-M 2 \ Sij5 k i, 



(A.17) 



n{gig 3 g k V u ) = M Q \ ijM + — - ( M 2 |y4<+ • • • ) + 



(5) 



M 4 | | OijdM+ 



d(d + 2) 



(2 A . n -1 

T 



-n- 1 [ M m M m + ••■ ) + ^— M 2 | ( M | ii <y w + — J + %rM 2 M [ 



(2) 

— m 2{0 \OijOkif ■■■ 



(9t9j9k9£) = 2{g i g j g k V u ) 
Summing over two repeated indices we get the two-index averages: 



i(9 2 VuVij) = n(V 2 gigj) = M 2 \ i3 + ~M 4 | £y + n 1 M 2 \ ^Mo\ij + ^-^2|0<%J 



(A.18) 
(A.19) 

(A.20) 
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n{g 2 9iVij) = n{(g-Vi)gigj} = M 2 \ij + -M^oSij +2n M \ ik M Q \ kj 



+— M 2 | 
a 



n(y?g % V l2 ) = n((g • Vi)ViiVy) - M 2{ij + -M Mo 5 tJ , 



1 n^ 1 
n((g • VOftVy) - M 2 |y + ^M 4 | <5y + n' 1 M 0{ih M \ kj + — M 2 , 



2M |y + ^M 2 | % 



Summing again, 



n(g 2 V?) =M m + n l M^ , 



"((g-Vi)^ 2 ) =M 4 | , 



n<(g • Vi) 2 ) = M 4 | + n^Mo^Moy, + -^-Af£ ]0 , 



n{(s ■ Vi).g 2 > = M 4 | + 2n- 1 M | ii M | J - j + n 1 ^^M^ , 



(. 9 4 )=2((g.V 1 )ff 2 ). 



A.21) 
A.22) 
A.23) 

A.24) 

A.25) 
A.26) 

A.27) 

A.28) 
A.29) 



Inserting Eqs. (| A.25|) - (| A.29|) into Eq. (|A.15|) one gets Eq. (f2T23|) . Analogously, from Eqs. (|A.14|) and (|A.13|) one 
jets, after some algebra, Eqs. (|2.24p and (|2.25p . respectively. 
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